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Abstract 

We study the geometric motion of sets in the plane derived from the homogenization 
of discrete ferromagnetic energies with weak inclusions. We show that the discrete 
sets are composed by a ‘bulky’ part and an external ‘mushy region’ composed 
only of weak inclusions. The relevant motion is that of the bulky part, which 
asymptotically obeys to a motion by crystalline mean curvature with a forcing 
term, due to the energetic contribution of the mushy layers, and pinning effects, 
due to discreteness. From an analytical standpoint it is interesting to note that 
the presence of the mushy layers imply only a weak and not strong convergence 
of the discrete motions, so that the convergence of the energies does not commute 
with the evolution. From a mechanical standpoint it is interesting to note the 
geometrical similarity of some phenomena in the cooling of binary melts. 


1 Introduction 

A definition of motion by curvature has been introduced by Almgren, Taylor and 
Wang (2] using a time-discrete approach as follows. Given a (smooth) set A 0 C 
as initial datum and a time scale r, one defines iteratively AJ. as a minimizer of 

A H > Per(A) -I— D(A,AJ._ 1 ), (1) 

r 

where AJ = Aq , Per(A) denotes the Euclidean perimeter of the set A and D(A, A') 
is a dissipation term that can be interpreted as the L 2 -norm of the distance between 
dA and dA'. The time-continuous piecewise-constant interpolations A T (t) = ATAj 
are then shown to converge to a time-continuous parameterized sets A(f), whose 
boundaries move by their mean curvature. Almgren and Taylor [T] have shown that 
the same scheme with a crystalline perimeter gives motion by crystalline curvature 
in dimension two. 
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The same scheme has been adapted to define a continuum motion for ferromag¬ 
netic spin energies on a square lattice whose static discrete-to-continuum approxi¬ 
mation is a crystalline perimeter by Braides, Gelli and Novaga [6J. In this process 
a scaling factor e > 0 and the corresponding perimeter Per e (A) for discrete sets 
A in el? are introduced, together with the corresponding discrete dissipations D e . 
These can be seen simply as the restriction of their continuum counterparts union 
A + [0,e] 2 of e-cubes. In this way discrete sets A £ k ’ T C el? are defined by iterated 
minimization of 

A H> Per 6 (A) H— D s {A,A £ k ff), (2) 

T 

where A £ f T are discrete interpolations of a continuum datum Aq together with 
the corresponding piecewise-constant-in-time interpolations A e ' T (t) = A?^ ,. The 
limit A of A e ’ T may depend on the mutual behaviour of e and r. 

This procedure can be framed in the theory of minimizing movements by De 
Giorgi (see [3]), which generalizes the approach of [2]. In [3] a notion of minimizing 
movement along a sequence of functionals has been given, that can be specialized 
for possibly inhomogeneous perimeter-type energies F e on el 2 : substituting Per £ 
with F e in the scheme above we can similarly define A e,T and obtain a limit motion 
A(t) passing to the limit both in e and r. In particular the following holds, upon 
the hypothesis of equi-coerciveness of F e and their T-convergence to some F: 

(i) ( pinning ) if r « e then A(t) = A 0 for all t > 0; 

(ii) ( commutation ) if e « r then A{t) is the minimizing movement of F with 
initial datum A 0 (hence, in the case of F e = Per £ the sets A move by crystalline 
curvature); 

(iii) ( critical scale) if e ~ r then the motion actually depends on e/r and is 
different both from (1) or (2). 

In [6] this last case is explicitly described through an effective motion: the limit 
Aft) depends on the ratio e/r, and this motion may depend on fine details of the 
energies F e and not only through their T-limit (see mi). The mechanism of 
evolution highlighted by the definition of A £ k ’ T is through local minimization of F e 
with a dissipation contribution that forces minimization on a small neighbourhood 
of the datum A k 'f 1 . If r is much smaller than e then by the discreteness of the 
parameters this neighbourhood contains the only A £ ff 1: and the motion is pinned. 
Conversely, if e is much smaller than r we can first pass in the limit as e —> 0 in 
([2]) and use the well-known property of convergence of minimum problems for T- 
convergence. In the critical case the motion optimizes the location of the interface 
combining energy and dissipation effects. 

In this paper we consider F e a sequence of double-porosity type energies, which 
mix ‘strong’ ferromagnetic interactions with ‘weak’ inclusions and still T-converge 
to a crystalline perimeter (see [5j). In this case the mechanism governing the time- 
discrete motion is of a different type from mm : since the perimeter energy due to 
weak inclusions in Aff, may be small with respect to the dissipation necessary to 
remove them from A k T , the latter is composed of a ‘bulky’ part, plus an external 
mushy layer composed of weak inclusions. This suggestive terminology is borrowed 
from theories in Fluid Mechanics where similar geometries appear in binary melts 
at solidification mmm- These mushy layer may then disappear at the next 
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step. As a result the final motion is not a simple motion by crystalline curvature, 
but it also contains a forcing term as a result of the effect of the mushy layers. 
This is true also for e « r for which we have the law of motion 

v = (an — b) + 

relating the velocity and the crystalline curvature. Note in particular that the 
conclusion (ii) above is violated. This is explained by a loss of coerciveness of the 
energies F e as e —>■ 0. 

2 Discrete setting and statement of the problem 

We are interested in describing a geometric continuum motion derived as the limit 
of time-discrete motions defined for discrete sets of eZ 2 as e —> 0 in the spirit 
of minimizing movements along a sequence of energies 0. In the specific two- 
dimensional case we are dealing with, the relevant information about the limit 
motion is obtained by considering initial data which are coordinate rectangles. 
The motion for more general sets can be derived from that case and obeys the 
same motion by crystalline curvature m E]- 

We will examine the time-discrete motions at fixed e and let eventually e —>■ 0. 
For the sake of simplicity of notation, we will state our problems in terms discrete 
subsets I in Z 2 without scaling them by e, but keep in mind that we are interested 
in the corresponding scaled sets el. 

With fixed a,(3 > 0, for e > 0 we define on I C Z 2 the energy 

F e (I) = e/3{(i,j) £ NN S : i £ /, j 0 1} + e 2 a {(i,j) £ NN™ : *£/,j^/}, (3) 

where 

NN S = {(*, j) £ Z 2 x Z 2 : ||* — jlloo = 1 and i\ = ji odd or * 2 = j 2 odd} 
NN^j = {(*, j) £ Z 2 x Z 2 : ||* — jlloo = 1 and i\ = j\ even or * 2 = j 2 even} . 

Besides these energies we introduce discrete dissipations, which take into account 
the L 2 distance of the boundaries of the discrete sets, given by 

De{I,I') = e 3 ( distoo(*,c(/ , ))+ ^ distoo(*,/')), (4) 

iei’\i iei\i’ 

where J, V C Z 2 and c(I') = Z 2 \ I'. 

Let e > 0 and r > 0 be fixed, together with a discrete coordinate rectangle 
Jo with the four vertices in (2Z 2 ) (we omit the possible dependence on e). We 
construct a sequence {/ ra }, where /„ minimizes 

E e (I,I n -i) =F e (/)+ 1 D e (I, In-i). (5) 

T 

The following lemma describes the relevant properties of the minimizers of (|5j) . 
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Lemma 1. For each n, if I n is a minimizer of ([5]) then 
In In— 1 ; 

b. there exist a discrete coordinate rectangle R n with the four vertices in (2Z) 2 
and a set W n C (2Z) 2 such that I n = R n U W n . 


This lemma (whose proof is postponed to the next section) shows that the 
discrete evolution of the sets If’ T = el n that we are interested in is described by 
a bulky part governed by the coordinate rectangles R e r f = eR n , and a mushy layer 
composed of weak islands given by Wf ,T = eW n . 

A relevant parameter in the description of the limit evolution is the ratio of 
time and space scales 


7 = 


lim 

£—>-0 £ 


( 6 ) 


where it is understood that r = t(e). 

We will show that the asymptotic description of the sets Wf ,T is not necessary 
to characterize the limit. Indeed 

• if 4ay < 1 then W n = ((2Z) 2 D/o) \ R n \ be., the mushy layer always contains 
all the weak sites in the initial datum Jo; 

• if 4ay > 1 then W n is contained in R n -i', he., the mushy layer at time step 
n—1 disappear at the next step. 

The case 4ay = 1 is exceptional, as in this case it may be equivalent in terms of 
the balance between energy and dissipation to maintain weak islands or ‘dissipate’ 
them. 

In any case, the relevant asymptotic description is given by the following result, 
whose proof is the content of the rest of the paper. In the description we do not 
treat in detail some non-uniqueness cases highlighted by the discontinuous right- 
hand side of the ODE in (101, which anyhow are completely analogous to those 
dealt with in El. 


Theorem 2. Let £ and r be fixed and let Rq be a given rectangle with the length 
of the horizontal side Lq and length of the vertical side Lq . Let Iq be the greatest 
discrete coordinate rectangle with the four vertices in (2Z 2 ) contained in fRn Let 
R n be the sequence of rectangles of Z 2 constructed by successive minimization as 
in Lemma [lj We define Lf and L 2 as the lengths of the horizontal and vertical 
sides of R n , respectively, and L\ T ft) = £L \ t / T , and L 2 £T (t) = eL 2 t ^ T j. Let e and 
t tend to 0 and (|6[) be satisfied; then L\ T (t) and L £T (t) tend to Lift) and L 2 ft), 
respectively, satisfying Li(0) = Lq and 


L[ (t) 


4 r 2/3y 
— maxi-— 

yL l 3 L 2 (t) 


2Q7 1 Pi 111 

3 6’2L 2 (*) 4 /J 


(7) 


for almost every t. In particular, we have pinning (no variation of L\) if 


respectively. 


4 B'y 

L 2 {t) > - -— z/4ay < 1, or 

4ay + 5 

L 2 (t)> 2 ^ < - if 4oy > l, 

The analog description holds for L 2 ■ 


( 8 ) 

( 9 ) 
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In terms of the crystalline curvat ure, which for a coordinate edge is given by 
k = (L being its length), equation (10 1 reads as 


2 

v — — 

7 L 


max' 


■ /3y k 
.~ 3 ~ 


2ay 

~Y~ 


1 /3jK 1 
6 ’ ~T~ + 4 


}J. 


( 10 ) 


where v is the velocity of the edge. This equation can be extended to coordinate 
polyrectangles and then to more general sets by approximation | 6 ]. 


Remark 3. Note that for 4aj < 1 equation (10) simplifies to 

2/37 


m = -~ 


3 L 2 {t) 


2 ay 


1 

6J 


( 11 ) 


Remark 4 (extreme cases). We can consider the cases r « e and e « r by 
letting 7 —» 0 and 7 — > + 00 , respectively. In the first case we have pinning for 
all initial data, and the motion is trivial. For 7 —» +00 the motion of L\ (and 
similarly that of L 2 ) is described by 


L[(t) 



U P \ 

2/3 1 

max | 7 

$U 2 (t) a ) 

’M*)J 


( 12 ) 


Remark 5 (a non-commutability phenomenon). In [5] it is shown that the discrete 
energies E e (I) = F e (-I) defined on subsets of el? T-converge to the crystalline 
energy 


F(A)= £|M|i dU\ 

JdA z 


whose minimizing movements give motion by crystalline mean curvature v = 

(see Almgren and Taylor 03 )> which corresponds in the case of a rectangle to side 
lengths satisfying 


L[(t) 


2/3 

L 2 {t) 


A general result in [3j shows that there exist a sufficiently slow time scale r such 
that the minimizing movement of a sequence of energies E e along r gives the 
minimizing movement of the T-limit. This is seemingly in contrast with the result 
in the theorem above since for 7 —» +00 (which corresponds to slow time scales) we 
have equation (12). This discrepancy is explained by the lack of equicoerciveness 
of the energies. Indeed the result in [4] only holds if the sequence is strongly 
equicoercive, and the appearance of the mushy region exactly corresponds to a 
weak (and not strong) convergence of the evolutions. 


3 Description of the structure of minimizers 

This section is devoted to the proof of Lemma [I] To show the result, it is useful to 
give a notion of connectedness for a discrete set / C 1? . We say that a discrete set 
/ C Z 2 is connected if the set a(I) = Uie/(* + [— 5 , |] 2 ) is connected. Given /, I' 
discrete sets with / C /', / is a connected component of I' if a(I') is a connected 
component of a(I). 
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Proof of Lemma [T} Note that setting R 0 = I 0 and Wq = 0 we have Iq = Rq U Wq. 
Suppose that such property holds for n — 1. 

If I n is a minimizer of |5]) and I n C (2Z) 2 , we can choose R n = 0 and W n = I n . 
If there exists i £ I n \ I n - i, then E e (I n \ i, i) = E e (I n , i) - 4e 2 a - e 3 r _1 < 
E e (I n , i). Since J„ is a minimizer, this implies /„ C J ra _! and the thesis follows. 

Now, we consider the case /„ \ (2Z ) 2 ^ 0 . Let C be a connected component 
of I n -1 such that C \ (2Z ) 2 7^ 0 . We show by contradiction that C C R n _ 1 . 
Since in a connected set I C Z 2 we have #(/ H (2Z) 2 ) < #(/ \ (2Z) 2 ) + 1, setting 
In = In\(C\ Rn-i) we get 

£ 3 

E s (I n , In- 1) ~ E e (I n , I n - 1) <-2/?£ 

T 

which is negative for £ small enough. Hence, C C R n -i- Denoting by R(C) the 
minimal discrete coordinate rectangle including C, since F e (R(C)) < F e (C), neces¬ 
sarily C = R(C) and C = ([a, a] x [ 6 , &])flZ 2 for some a, a, 6 , b £ Z. We show that the 
four vertices of C belong to (2Z) 2 . Reasoning by contradiction, it is not restrictive 
to assume C = ([a, a] x \b, 6 ])flZ 2 with a ^ 21. Note that, setting N = #([&, b] nZ), 

we have Yn=b d (( a ’ *)> c ( J n-0) > 2 -Ne - 2 - Since rf (( a > *)> c(/„_i)) > 2/3 t£ - 2 , 
setting I n = I n U ({a — 1} x ([ 6 , 6 ] D Z)) it follows that 

3 ^ 

E e (I n ,I n -i) - E e (I n ,I n - 1) < 2«£ 2 --V d((a-l,i),c(I n _i)) 

T z —' 

i—b 

< 2aff2 “ ^ ' 52 d((a,i),c(I n -i)) + £ — 

i=b 
£ 3 

< 2a£ 2 — /?£ H- 

T 

which is strictly negative for £ small enough. 

Now, we show that the connected component C = ([a, a] x [fo, 6]) fl Z 2 is the 
unique connected component intersecting Z 2 \(2Z) 2 . We prove this by showing that 
the center ( 711 , 712 ) of R n -i belongs to C. Since C \ (fill) 1 7 ^ 0, it is not restrictive 
to assume b <b. If a < m, we consider the set I n = I n U ({d + 1} x ([&, 6]) n Z). 
Since 

b 

£ 3 t -1 ^2 d((a, i),c(I„.)) > 2 e/3 + 6£ 2 a, 

i—b 

we get 


b 

E e (I n , I n -i) - E e (I n , J„_i) < -£ 3 T -1 d((d, 7 ), c(/ n )) + 2 e/3 

i—b 

< -12£ 2 a < 0. 

The same argument holds for a > n\. Hence, a < n\ <a. Assume by contradiction 
that b < ?r 2 , and define I n = I n U (([a, a] fl Z) x ([6 + l,n 2 ] H Z)) (if n 2 fL 2Z, we 
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substitute 712 by 712 + 1). Since 

3 ^ 

{d{(i,b), C (I„)) + d((i,b - 1 ),c(I n ))) > 2s/3 + 2e 2 a, 

i=a 

we get, recalling that d((m, n 2 ), c(7„)) < d{[n\, b), c(7„)), 

3 ^ 

E e (I ni I n - 1 ) - E e {I n ,I n - 1 ) < -A'^- ^ (d((i,7),c(7„)) + d((i,b - l),c(7 n ))) 

i—a 

+2Ke/3 + 2 Ke 2 a 

< 0 

for some K > 0. Hence, b > 712 . The same argument shows that 6 < 712 , and the 
claim is proved. Since I n cannot contain isolated points in (2Z) 2 which are not in 
In— 1 , the proof is complete. □ 


4 The iteration procedure 


Given L, L' > 0 we define 

i e = ([o,* e (L)] x [o,i £ (i')])nw 2 

where i £ (:r) denotes for any x > 0 the greater even integer less than |_§ J; that is, 


i e (x) = 2 


b/gj 

2 


J- 


Setting 7° = 7 e , for 71 > 1 we denote by 7" a minimum point for the energy 

A e (7,7" _1 ) = F £ (7) + ^-D e (I, 7" _1 ). 

Setting for (7,7) £ [0,i £ (Z/)/4] x [0,7 £ (L)/4] 

C £ (7,7) = [27,i e (L') - 27] x [2 k,i e (L) - 2k], 


(13) 


(14) 


(15) 


Lemma [I] ensures that if 7 is a minimizer of (14) then there exist (7", 7") £ 
([0,i e (Z/)/4] x [0, z e (L)/4]) D N 2 and C N 2 with C e {K,k?) n N 2 C W” C I e 
such that 


7 = 7 e (7", 7", W”) = (C £ (7 £ , 7") n N 2 ) U (IT” n (2N) 2 ). 

Moreover, 7£ > 7™" 1 , 7£ > 7£ t " 1 and W e n C W^" 1 . 

In the following sections, we show that, up to a subset of (2N) 2 , the minimum 
problem for the energy E e (I,I e ) has a unique solution 

II = ( C e {h\,kl) n N 2 ) U (Wl n (2N) 2 ) 


7 





with h\ and k\ independent of e for e small enough. Hence, the corresponding 
result holds for In particular, we show that if the initial length L of an edge is 
lower than a critical threshold then the corresponding “displacement” h depends 
on the initial length as follows: 


h = h(L) 


2/?7 2ay 1 | 
-~3lT ~~ + 6-T 


As to the set Wf of the weak islands (i.e., points with weak connections), note 
that the variation of the energy due to an isolated set {i} C (2N) 2 is given by 


4ae 2 — —distoo (i,C(/ e )). 

T 


Hence we may rule out the formation of weak islands only if 4a — - > 0, that is if 


4ay > 1. 


(16) 


Thus, we consider two different cases in dependence of the value of ay. 


5 Case 4cry < 1 


Recalling the properties shown in the previous section, it follows that in the mini¬ 
mum problem for the energy E e (I,I e ) defined in ( |14[ ) we can consider only sets of 
the form 

I e (h, k ) = (C E (h, k) n N 2 ) U (J E n (2N) 2 ) (17) 

for {h, k ) € ([0, x [0, ^]) n N 2 , and the minimum problem for E e corresponds 
to minimize 

fe(h,k) = 1 (F e (I e (h,k)) - F e (I e ) + -D s (I e (h,k),I e )). 

Since we are interested in e —> 0, it is not restrictive to consider the case when 


min{L, L'} 
2 n 


for n £ N large enough 


(18) 


(hence i e (L ) = j if L < L' or i e {L') = if L > L'). We will then make this 
assumption, commenting on the error that we make under this hypothesis when 
necessary. 

In the sequel of the section, we prove the following result. 

Theorem 6. For e small enough, the energy E e (I e (h, k), I e ) has a unique mini¬ 
mum point in ([0, |^] x [0, jr]) nN 2 given by 


( 0 , 0 ) 

{n(L),n(L’)) 


if min{T, L'} = X c < max{A, L 1 } 
otherwise 




where for any l > 0 


»(0 = 


0 if l > \ c 

L^-^ + gJ 

and the critical threshold X c is given by 

x _ 4/37 

Ag — . 

4a7 + 5 

Proof. Suppose L' > L. First, we note that, in the case L' > L , if 4 eh > L then 
the variation of the energy is strictly positive (for e small enough), so that the 
minimum point does not belong to this set and 


(19) 


( 20 ) 


fs(h,k)> min f e 
[o,4] 2 nN 2 


if h> 

4e 


( 21 ) 


Indeed, the dissipation term turns out to be larger than 

a 3 /IA/( 2e )J / j t 

* ( £ (E'+A sH 

\ j=i N=i 

Li/(4e)J / j T 

- £ (£(2l-l) + (2i-l)([ 5 J-j) 

j =i v ;=i 

and thus, up to a uniformly bounded term, larger than . Since the variation 
of the boundary energy F e is uniformly bounded, for e sufficiently small it follows 
that 

f e (h, k ) > e~ 2 C if L < 4eh < L' 


where C > 0 is independent of h,k,e, and (21) follows. 

Now, we consider (h,k) £ [0, ^] 2 ON 2 . Recalling that by assumption £ 2N, 
the variation of the boundary energy is given by 


F e {I e {h 1 k)) — F e (I e ) = 4e(— [3h + aLh — ft k + aL'k) 

+4 ae 2 (h + k — 4 hk) — 4 ae 2 g e k 


( 22 ) 


where g e = h — i e (I/) £ [0,2). Note that if we do not assume (18) then we have 
an additional term —4 ae 2 Q^h with ^ — * e (L) 

As for the dissipation term, setting M = max{/i, fc} we have 


D(I e (h,k),I e ) = 2 £ 3 ((- + 1 - 4Af) £ j - (- + 1 - 2M) £(2j - 1)) 

3=1 3=1 

+2 £ 3 ((z £ (L') + 1 - 4M) + 1 - 2M) ]T(2j - 1)) 

3=1 3=1 

2M 3 M 3 

+ 4e3 (E (E 1 + i( 2M - ■?')) - E (-!) + ( 2 -? - !)( M - j 


3=1 1=1 


3=1 1=1 
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Hence, 


-D(IJh,k),I e ) = -(-h(3h + 2) + 2h 2 + 2h-12h 2 M-8hM) 

T T V £ / 


c-3 / tJ v 

+ — ( — fc(3fc + 2) + 2fc 2 + 2fc - 12k 2 M - 8 /cMj 

3 3 

- — (3 g E k 2 + 2 g e k) + — ( 8 M 3 + 8 M 2 ) 

T T 

= -(Lh(3h + 2) + L'k(3k + 2)) 

e 2 

+ —( 2 h 2 + 2 h + 2 k 2 + 2k - 8 hk - 12Mmm{h, k} 2 - 4M 3 ) 
7 

e 2 

- g e k(3k + 2). 

7 


Again, if we do not assume (18) we obtain an additional term — — g^h(3h + 2) with 

Qe = j- *e( L )- 
Setting 


7r(x) = 3ar + 2(2(07 + l)x 


and, for l > 0 , 


Pi(x ) = —tt(x) — 4 fix = —x 2 — 2(2/3 — (2a H—V) 3 ' 
'Y O' V V 7/ / 


(23) 

(24) 


we get 


/ e (h, fc) = Pi(h) + P L ,(fc) + -P(/i, fc) - -gMk) 

7 7 

where R{h , fc) is the symmetric function defined for h > k by 

P(/i, fc) = 2 /i 2 + 2 ( 2 <o 7 + l)/i + 2k 2 + 2(2 a 7 + l)fc 

— 8(207 + l)hfc — 12 hk 2 — 4 ft 3 . 


(25) 


Without assuming (18 1 we get an additional term — ~g^Tr(h) in (251. 

The sequence {/ £ } uniformly converges on the compact sets of R 2 to / defined 
by f{x,y ) = Pl{x) + Pu{y)- We denote by (m(L),m(L')) the minimum point of 
/ in R 2 , so that m(l) is defined by 

2/3y — (2a7 + 1)£ 


»(0 = 


3/ 


(26) 


Note that if L = P' the function / e is symmetric. If L' > L, then the symmetry 
of R gives f E (x, y) - f e (y, x) = (P L - P L - - - (Pl - Pl> - y-7r)(y); since 

Pl — Pl' — — 7 T is strictly decreasing in ( 0 , + 00 ) for e small enough, it follows that 


fe{x,y)> f e {y,x) if 0 < x < y. (27) 

Hence the minimum of f e in [0, jz } 2 is achieved in {(x,y) : x > y}. 

Now we prove that the minimum of f e is in fact achieved in a compact set 
independent of e. 
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Proposition 7. There exists x > 0 independent of e such that, fore small enough, 
if (x E ,y e ) is a minimum point of f e in [ 0 , jr] 2 , then ( x e ,y E ) £ [0,5i] 2 . 

Proof. Since f e coincides in {a; > y} with a polynomial function of degree 3, then if 
L < L' (hence m(L ) > m{L')) the uniform convergence ensures that the (unique) 
critical minimum point of f e belongs to a compact neighborhood of (m(L),m(L')) 
included in {x > y}, and independent of e. In the case L = L' , the computation 
of the partial derivatives gives V/ £ ^ (0,0) in [0, j }] 2 n {x > y}. 

We fix x > ma x{m(L) + 1,0}; the minimum of f e in [0, jf] 2 n {x > x, x > y} 
is then achieved on the boundary, where a computation shows that 


fe(x,y) > Pl{x) + min P L , + o(l) e _> 0 . 


Choosing x such that PiXx)+\mn Pjj > 0, then in [0, jf ] 2 n{a" > x, x > y} we have 
f e (x, y) > 0 = / e (0,0). Since f e (x, y) > f e (y, x) if x < y, the thesis follows. □ 

To prove Theorem [oj we hav e to consider the minimum problem for f e in 
([0,£] x [0, ^]) fl N 2 . Recalling (21), the minimum is achieved in [0, |] 2 fl N 2 


and thanks to Proposition [7J it is sufficient to show that f e has a unique minimum 
point independent of e in [ 0 ,x ] 2 D N 2 . 

If / has a unique minimum point in [0,5 ;] 2 PlN 2 , then the uniform convergence 
ensures that, for e small enough, it coincides with the (unique) minimum point of 
f e in [0,5 ;] 2 fl N 2 , concluding the proof of the uniqueness and independence on e. 

The minimum of / is achieved in the points (h, k) in [0, ir] 2 HN 2 which minimize 
the distance from (m(L), m(L')). Thus, the uniqueness fails only if m(L) + } or 
m(L') + \ belongs to N\ {0}. In this case, we have to compare the values of f e — f 
in the set of the minimum points of /. 

Note that the condition m{L) + \ < 1 corresponds to L > A c , where A c is the 
critical length defined in (20 1 . Hence, recalling that for l > 0 


m 1 2 /37 2c*7 1 

mW + 2 = Yi —r + 6- 

if both m{L) + A and m(L') + A do not belong to N \ {0} then the minimum point 
of f £ in [0,S] 2 nN 2 is unique and it is given by 

- (0,0) = (n(L),n(L')) if L > A c ; 

- ([ m{L) + \\, 0) = ( n(L ), n(L')) if L < A c < L'\ 

- ([ m{L) + Aj, [ m(L') + Aj) = (n(L),n(L')) if V < A c 


where the function n is defined in (19). 


It remains to check the case when m{L) + A or m(L') + A belongs to N \ {0}. 
Note that R(x,y ) — g e ^(y) decreases with respect to each variable in {x,y > 1}, 
and also R(x, 0) is strictly decreasing in {s > 1}; moreover, for any x > 1 then 
R(x , 0) > R(x, 1) —p e 7 r(l). The monotonicity properties of R(x, y) — g e r’(y) allow to 
deduce that the minimum is achieved in (| m(L) + , [m(L') + |J) = ( n(L ), n(L')) 

for each L,L' except when L = A c . In this case, if also L' = A c , the comparison 
between / e ( 0 , 0 ), / e (l, 0 ), / e ( 0 , 1 ) and / e (l,l) ensures that the minimum point is 


11 





(1,1), again corresponding to ([m(L) + |J, [m(L') + lj). If L = X c and L' > A c , 
since f e ( 1,0) = 4ea > f E ( 0,0), then the minimum is obtained in (0,0) = ([m(L) — 
WAm{L') + \\). 

This completes the proof of Theorem [6j □ 

In particular, except for the case when one or both the initial lengths are greater 
than A c , the minimizing set for the energy E e is 


I e (n(L),n(L')) = J e ( 


2/37 2«7 1 
”3Z 3 h 6J 


2/37 2«7 1 


L3L' 


6 


J) 


with I E defined in 

We can use the set I e (n(L),n(L')) as a recursive datum for the energy (14). 


Note that the weak sites of the initial configurations will be part of each minimizers. 
We can then describe the motion through the velocity of the moving rectangle 
corresponding to C e (h, k). 

We finally compute the velocity of a side. Consider an edge with initial non- 
scaled length L. We proved that it moves only if L < A c (also for L = X c if L is 
not the minimal length of the edges). In this case, the displacement of the edge is 
given by 

Wl 2«7 . 1 I (28) 


3 L 


sJ- 


Hence, if L < X r 


v 7 (L) 


2 2/3y 2a7 1 

7 L"3L" ~ “3“ + 6- 


where r(L) = m(L) + \ — [m(L) + |J. 


4/3 t 4a 6 r(L) - 1 

3 L + Y + 3t 


6 Case 4ct7 > 1 


The properties of the energy shown in Section [4] imply that a weak island may 
appear only when the distance from the boundary is strictly greater than 4ay. 
Hence, except if 4aj is an odd integer, an isolated point i G (2N) 2 belongs to a 
minimizing set for E e (I , I e ) if and only if the distance from the boundary of I E is 
greater than the minimal odd integer larger than [4a7j, which is given by 2N aj — 1 
where 


N — 

1 y Q7 — 


|_4a7j + 1 
2 


Recalling the definition (15) of C e (s,t), if we define for (s,t) £ [0 ,i s (L')/4] x 
[0, i E (L)/4] 


J E {s,t) = (c e (M)nN) U (c e (A7 T ,A7 7 )n(2N) 2 ) (29) 


where we omit the dependence on L and L ', then if 4a7 is an odd integer it follows 
that 


E e (Je(s,t), I E ) = E e ({C e (s,t) D N) U (W n (2N) 2 ),/ e ) 


(30) 
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where W C C e (N ai — l,iV a7 — 1). Hence, except if 4cry is an odd integer, the 
minimum problem for the energy E e (I,I e ) corresponds to minimize 


9TM = ~(E e (J e (s,t),I s ) -E e (I e )) 


(31) 


in ([0, i e (L')/ 4] x [0, * £ (L)/4]) n N 2 . 

In the following theorem we give an explicit computation of the set of minimum 
points of .g“ 7 (s,t), showing that this set is independent of e (for e small enough), 
and that there is uniqueness except if 4 q?7 = 2 and min{ L, L 1 } = A + < max{T, L'}. 
We prove the following result, which takes care of almost all values of 4ay, except 
when this is an odd integer. In that case we may have the non uniqueness phe¬ 


nomenon as described in (30), which will not affect the final description of the 
motion. 

Theorem 8. Let 4aq £ ffi. \ (2Z + 1) (i.e., not an odd integer). If 4aq > 2, then 
the unique minimum point of g < f 1 {s,t) in the set ([0, i e (L')/A] x [0, z e (L)/4]) D N 2 
is (tp(L),tp(L')) where 


<P(l) = 


and the critical thresholds are given by 

4/37 


0 

if l> X+ 

\ftL + 1 I 
l 21 ' 4 J 

if l e (X~ t A+] 

3Va 7 

ifl€(X*,X~] 

| 2/37 2<27 , 1 1 

L 31 3 ' 6 J 

if i<K 


A! = 


4a7 + 5 + 6 N, 


X~ = 


2/3 7 


O'7 


-1 


and 


x + = 2 -^. 


(32) 


(33) 


If 4«7 < 2, the minimum point is again unique, and it is given by 
j (0,0) if min{L, L'} = X + < max{L, L'} 

I ((p(L),<p(L r )) otherwise. 


If 4«7 = 2, for min{L, L 1 } = A + < max{L, L'} there is no uniqueness, and the 
minimum points are (0,0) and = (1,0). Otherwise, the minimum 

point is unique and it is given by (ip(L),ip(L r )). 


Remark 9. (i) Note that the non-uniqueness exactly for 4ay = 2 depends on the 
simplifying choice (18). This choice however influences only the value at which we 
have non uniqueness, but not the final description. 


(ii) Note that ip can be equivalently written as the integer part in (10). 


Proof. We assume L' > L. Since the weak islands appear only when the distance 
from the boundary is greater than 2 N aj , in the computation of the minimum point 
of the energy we have to consider three cases in dependence on the values of s and 
f, namely s,t < N a7 , min{s,f} < IV Q7 < max{s,f} and s,t> N ai . 
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We start by computing the expression of g“ 7 in [0, N ai ] 2 n N 2 . 

Assuming s > t, for (s,t) € [0, N ai ] 2 flN 2 we get 

4 2s J 

g“ 7 (s, t) = —4 ps — Apt — 4 eas — Aeat + — EE /+j(2s-j)) 

7 j=i i=i 

9r 2s p F 2 * 

+ -(i e (L) + 1 - 4s) £ j + -(* e (L') + 1 - 4s) £ j 

7 j=i 7 j=i 

L U e 

= —p(s) — ApS H- p(t) - Apt --g e p{t) 

111 

+ ~(^(! + 4s)p(s) + (1 - 4s)(p(s) +p(t))) - 4a 7 (s + t)j, 

where p(x) = 2x(2x + 1) and g e = — i e (L'). Hence, setting for l > 0 

Qi(x) = —p(x) — APx = — x 2 — 2^2 P — ~^jx, (34) 

the function g“ 7 can be expressed in [0, N ai ] 2 n N 2 as 

5e 7 (M) = Ql(s) + Q L '(t) + —r(s, t) - -g £ p(t) (35) 

7 7 

where r(s,t) is the symmetric function defined for s > t by 
4 

r(s,t) = -(1 + 4s)p(s) + (1 - 4 s)(p(s) +p(t)) - 4ay(s + t). 

The following lemma holds. 
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Lemma 10. //4«7 > 2, the unique minimum point of g < f 1 {s, t) in the set [0, iV Q7 ] 2 n 
N 2 is given by (<p(L) A N aj , f(L') A iV a7 ). 

If 4a7 < 2, the minimum point is again unique, and it is given by 
I (0,0) if L = A+ < L' 

( (<p{L) A N ai , ip(L') A N ai ) otherwise. 


If 4a7 = 2, for L = X + < L 1 there is no uniqueness, and the minimum points are 
(0,0) and (<p(L) A N ai ,ip(L') A N ai ) = (1,0). Otherwise, the minimum point is 
unique and it is given by (ip(L) A N a ^,ip(L') A JV a7 ). 

Proof. The result follows by the uniform convergence of er(s, t) — eg e p(t ) to 0 in 
[0, -/Vq, 7 ] 2 , so that the minimum of the energy is obtained in the set of the points in 
[0, N ai ] 2 flN 2 minimizing the distance from the minimum point of Ql(s) + Quit), 
given by p.(L')), where 

2/3y - l 

M0 = —^—■ (36) 

This set contains only one point except if p(L) + \ or p{L') + \ belong to 
the set of integers {1,..., N ai — 1}. In all other cases, noting that the conditions 
p(l) = N aj — \ and p(l) = \ give the critical lengths A - and A + respectively, the 
thesis follows. 

When the set of minimum points of Ql(s) + Quit) in [0, N ai ] 2 n N 2 contains 
more than one point, we need some monotonicity properties of r. The function 
r(s,t) — g e p(t) is strictly decreasing with respect to each variable in {s,t > 1}; 
moreover, in {s > 1} the function r(s, 0) is strictly decreasing and r(s, 0) > r(s, 1) — 
g e p( 1). This gives the thesis except for the case p(L) = \ and p.(L’) < \ (hence 
L = X + and L' > A + ). By noting that 

r(l, 1) - g e p{ 1) < min{r(l, 0); r(0,1) - g e p( 1); r(0,0)}, r( 1,0) = 2 - 4^7 


the proof is complete 


□ 


Now, we compute the expression of g“ 7 in ([iV a7 , |^] x [iV a7 , ^]) n N 2 . 

Following the proof of Theorem [6j it turns out that it is not restrictive to 
consider only (s, t ) € [iV Q7 , ^] 2 n N 2 . 

We can decompose </“ 7 (s,t) in the sum of two contributions. The first is due 


to the set J e (N c 

A* 7 ( N .'< 


j N ai 


9s 


a7 , iV a7 


) and it is given by (35) 
) = Qi(3V a7 ) + Ql>{N a 7 


. £ 

) + -r 
7 


(A ^ a7 , A (q , 7 


) QeP{N a j) ■ 

7 


Then, we have a contribution which can be computed as in the case 4ay < 1 
by substituting the initial set I e with J“ 7 = C e {N ai , N ai ) n N 2 . Hence, following 
the proof of Theorem [6j the contribution to gf 1 is given by 

/“ 7 (s - N ai ,t- iV Q7 ) + - 2 N ai #(/“ 7 \ J e (s,t)), 

7 
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Figure 2: The set J £ (s,t ) with s,t>N ( 


aj 


where the term 2iV Q7 #(/“ * * * * * 7 \ J E (s,t)) takes into account the additional distance 

2 N ai of each point of the set from C(I E ), and 


1 


k) = - [E e {J e (h + N ai , k + N ai ), Ir) - F e (ir) 

= -P L _ 4eA r (/i) + PL'-ieN^{k) + —R(h, k) - -g e Tr(k), 

7 7 


(37) 


with 7r, Pi and R defined as in (23), (24) and (25) respectively. Since 


P L -ieN a ,(h) = P L {h) --(12N ai h 2 + 8N ai (2a^ + l)h), 


we get for g“ 7 (s,f) the following expression 


9 r(s,t) = P L {s — N ai ) + 


6 LN, 


a. 7 


(s — fV a7 ) + QL{N ai ) 


7 

+-Fz/(i — N a ~f) -\ -— (t — N<x~i) + Qu{E[ ai ) 

7 

£ £ £ 

~\ r(lV a7) TV a7 ) -)- R a j(s N a <y,t N af y^ 

7 7 7 

where 7r a7 (fc) = 3A 2 + 2(3_/V a7 + 2ay + l)fc + p(lV a7 ) and 

R ai {h, k ) = i?(/i, A:) + 41V 0 , 7 (1 — 6 N ai )(h + k) — 24 N ai hk 
— 12N aj (h 2 + k 2 ) - 8N aj (2a'y + l)(h + k ). 

Now we prove the following lemma. 

Lemma 11. The function g £ ' y (s, t) has a unique minimum point in the set ([iV a7 , |^] x 
[ N on, ^DnN 2 , given by (tp(L) V N a ^,ip(L') V iV a7 ). 
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Remark 12. Note that ip(l) V N ai = [^§r~ + \ J for A < A* and ip(l) V N ai = 


N ai otherwise, where A* is defined by (33). 


Proof of Lemma 0 We note that the minimum point of Pi (a:) - 


6 lN a 


- [x) is given 


by m(l) — N ai , where m(l) = — pp — A denotes the minimum point of P;; the 

condition m(l) — N a7 = \ introduces the critical threshold A* given by (33). 

As in Proposition [7] we can prove that also in this case the minimum of gf 1 
is in fact achieved in a compact set independent of e. This allows to show that a 
minimum point of in [N aj , j^] 2 n N 2 necessarily belongs to the set of points 
minimizing the distance from (m(L),771(1/)). 

Since the function R a ^(x,y) — g e Tr aj (y) is strictly decreasing with respect to 
each variable in {x,y > 0}, then the thesis of Lemma [TT| follows also if the set 
of points in [fV a7 , -p] 2 n N 2 minimizing the distance from (■ m(L),m(L ')) contains 
more than one point. Note that the minimum is obtained for (\m(L) + Aj,0) even 
if L = A* and L' > A*. □ 


Now, we have to consider the case min{s, t} < N a ~ < max{s, t}. 

If (s,t) € ([lV a7 , x [0, iVo, 7 ]) nN 2 , we can decompose g“ 7 (s,t) in the sum of 
three contribution. The first is due to the set J e (N ai ,t) and it is given by (35) 


g^ 1 (N ai ,t) = Q L {N aj ) + Ql,'( t) + -r(N aj ,t) - -g s p(t). 

7 7 

The second term can be computed following the same argument of the case s, t > 
N aj , and it is given by 

/“ 7 (s - iV Q7 ,0) + £ - 2N ai (3(s - N ai )i e {L) + 2(s - IV Q7 )(1 - 6 N ai )) 


where /“ 7 is defined as in (37). Moreover, we have to consider the dissipation 
due to J s (N aj , t ) \ J e (s, t); denoting by Q e (s, t) the set [2TV a7 , 2s) x [2 1, 2 N ai ), the 
additional term turns out to be 


2N a -,~2t 

- Y, dist °°(b C &)) = -(2s - 2A r Q , 7 ) £ (j + 2 1) 

^ Q e (s,t)ni e ~ i=i 

= —(s N ai )(N ai — t){2N ai + 2 1-\- 1). 
7 


Hence, we get for the energy in {[N al , p] x [0, A a7 ])nN 2 the following expres¬ 


sion 


g“ 7 (s,t) = P L (s-N ai ) + 


6LN a 


-(s — N aj ) + QL,{N ai ) + Q L> {t) 


£ £ 

+ -r„ 7 - - g E p(t ) 
7 7 


17 








L 



a ® 

■ ffl 

■ H 

■ ® 
a a 
a a 
a a 


2N (ty 2s 


L' 


Figure 3: The set J £ (s,t ) with s > N ai > t 


where Pl is defined by (24) as in the case 4«7 < 1, and 


r a ~f{h 1 t ) — r(N ai , t) + 8h(N a j — t)(2N a j + 2f + 1) + -R Q7 (/i, 0). 


Note that r ai (h,t) decreases with respect to each variable in {h,t > 0}. 

The following lemma holds. 

Lemma 13. If L < A*, the unique minimum point of g’f 1 (s,t) in the set ([IV Q7 , 

[ 0 , JV a7 ]) n N 2 zs (^ + 1\ ML') A iV a7 ). 

The proof follows as in the previous cases by uniform convergence; the minimum 
is obtained in the set of the points in ([lV a7 , L.\ x [0, lV a7 ]) n N 2 minimizing the 
distance from (m(L), /i{L')). 

Again, the monotonicity of r a7 implies that the result holds also when m(L) + 1 
or p(L') + \ are integer. 


If ( s,t ) € ([0,iV Q7 ] x [N ai , T]) n N 2 , following the argument of the previous 
case the expression for g“ 7 turns out to be 


9 r(s,t) = Q L (s) + P L '(t-N a7 ) + ^^(t-N ai ) + Q L ,(N ai ) 

7 

£ £ 

H r a7 (t N a7 , s) Pe 2 r a7 (t) 

7 7 

where 7r Q7 (fc) = 3 k 2 + 2(3_/V a7 + 2ay + 1 )k + p{N ai ). 

The following lemma holds. 

Lemma 14. If L' < A*, the unique minimum point of g^' y (s,t) in ([0,-ZV a 

[N* 7 , if}) nN 2 (N aj , L|^- ? f r + gJ) = ML) A N a „ - *p. + ij). 
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Remark 15. If L > A*, then the minimum of g“ 7 (s, t ) in ([iV Q7 , x [0, IV a7 ])nN 2 
is achieved for s = N al . This implies that the minimum point of the energy in 
([0, x [0, IV a7 ]) nN 2 is unique and belongs to [0, N ai ) 2 nN 2 . The corresponding 
result holds if L' > A*, so that in this case the minimum point in ([0,A a7 ] x 
[0, j:]) fl N 2 is unique and belongs to [0, iV a7 ] 2 n N 2 . 


The thesis of Theorem [ 8 ] follows by comparing the results of Lemmas 10 13 
m an d and recalling Remark [15] Q 


We can use the set I e (n(L),n(L')) as a recursive datum for the energy (14). 
Note that the weak sites of the initial configurations not in R n will disappear at 
the following step if 4aq > 1 , while the case 4aq = 1 is exceptional and we may 
keep or discard any of such weak sites in the following. In any case we can describe 
the motion through the velocity of the moving rectangle corresponding to C e {h , k). 

We finally compute the velocity of a side. We consider an edge with initial 
non-scaled length L. Theorem [ 8 ] shows that it moves only if L < A + (also for 
L = A + if L is not the minimal length of the edges, or if 4ay > 2). In this case, 
the displacement of the edge is given by 2 <p(L) Hence, the velocity of each side is 
given by 

2 


f (L) = — y(L). 


Remark 16 (the limit case 7 


Too). Since lim 

7—>00 


2N a . 

7 


a: = a:( 7 ) = 


4/3 7 


= 4a, we have 


P 


4a 7 T 6 N ai + 5 4a 


Choosing 

we get the limit for 7 —> Too 


L < 


4a 


4 8 4 a 

lim v L ( 7 ) = --T-. 

7—>+00 6 L o 
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